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We analyse available experimental data on the total charged-current 
vN and VN cross sections for quasielastic scattering and single-pion neu- 
trinoproduction. Published results from the relevant experiments at ANL, 
BNL, FNAL, CERN, and IHEP are included dating from the end of six- 
ties to the present day, covering and beams on a variety of nuclear 
targets, with energies from the thresholds to about 350 GeV. The data are 
used to adjust the poorly known values of the axial masses. 

PACS numbers: 12.15.Ji, 13.15.+g, 14.20.Gk, 23.40.Bw, 25.30.Pt 

1. Introduction 

It is well known that the theoretical description of the cross sections for 
CC and NC (quasi)elastic neutrino-nucleon scattering (QES) and single-pion 
neutrinoproduction through baryon resonances (RES) are very sensitive to 
the shape of the weak axial-vector elastic and transition form factors. By 
adopting the standard dipole parametrization for these form factors, their 

OPS 

shapes can be described with the two phenomenological parameters 
and , the so-called axial (dipole) masses. In general, these masses are 

different and, moreover, the numerical value of is vastly dependent 

of the particular dynamic model for the resonance production. 

The experimental values for both M^'^ and coming from mea- 

surements of (quasi) elastic neutrino and antineutrino scattering off pro- 
tons and nuclei and from the more involved and model-dependent analyses 
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of charged pion electroproduction off protons, spread within rather wide 
ranges. In this study we attempt to fine-tune the axial masses by fitting all 
available data on the CC QES (with Ay = 0) and RES Ivr total cross sec- 
tions for and scattering off different nuclear targets from experiments 



at ANL [1-5], BNL [6,7], FNAL [8-12] CERN [13-26], and IHEP [27-33]. 



In our opinion, this procedure is more selfconsistent in comparison with the 
usual straightforward averaging over the experimental values of the axial 
mass (see, e.g., Ref. [34]) extracted under different assumptions about the 
other badly known factors involved into the analyses of each experiment. 

2. Axial mass from the data on quasielastic scattering 

Figure ^ shows a compilation of the QES data from experiments at 
ANL [1,3,4], BNL [6], FNAL [10-12], CERN [13-17,20,26], and IHEP 
[27-30, 32, 33] performed with a variety of nuclear targets. The cross sec- 
tions reported in the earlier experiments [1,3,13-16] exhibit uncontrollable 
systematic errors and fall well outside the most probable range determined 
through the fit to the full dataset of about 200 datapoints; the value of 
evaluated for each subset of these data exceeds ~ 5 ndf. Hence, following 
the (nonstringent) selection criterion x^/ndf < 4.5, they were excluded from 
the final fit. 

For the i^n — > fi^p and Vp fx^n cross sections we use the result 
of Ref. [35] neglecting possible second-class current contributions (see Ap- 
pendix); under this standard assumption it coincides with that of Strumia 
and Vissani [36]. For the elastic electromagnetic form factors we apply the 
QCD VM model of Gari and Kriiempelmann [37] extended and fine-tuned 
by Lomon [38] ("GKex(02S)" version) and the most current inverse polyno- 
mial parametrization by Budd et al. [39] ("BBBA2006") obtained through 
a global fit to the world data on the Sachs form factors. For the axial and 
pseudoscalar form factors we use the conventional representations [40] 



with Fa{0) = gA = -1.2695 ± 0.0029 [41] (assuming gv = 1) and Ma = 
j^QES ^jgjj^g ^ £j.gg parameter of our fit. 

The nuclear effects for the data obtained for deuterium [3, 4, 6, 10, 26] 
and neon- hydrogen [11] targets were subtracted by the authors of the ex- 
periments. Therefore these data are fitted by the cross sections evaluated 
for free nucleons. To describe the remaining experimental data we apply the 
relativistic Fermi gas model by Smith and Moniz [43] with the kinematics 
and values of binding energies and Fermi momenta of the target proton and 
neutron determined by the composition of each target quoted in Fig. ^ 
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Fig. 1. Total quasielastic v^n and Vf^p cross sections measured for different nuclear 
targets by the experiments ANL 1969 [1], ANL 1973 [3], ANL 1977 [4], BNL 1981 
[6], FNAL 1983 [10], FNAL 1984 [11], NuTeV 2004 [12], CERN BC 1965 [13], CERN 
HLBC 1966 [14], CERN HLBC 1967 [15], GGM 1973 (Gargamelle, CERN) [16], 
GGM 1977 [17], GGM 1979 [20], BEBC 1990 (CERN) [26], IHEP SKAT 1981 [27], 
IHEP 1982 [28], IHEP-ITEP 1985 [29], IHEP SKAT 1988 [30], IHEP SKAT 1990 
[32], and IHEP SKAT 1992 [33]. The curves and bands correspond to the world 
average value of M^'^ = 0.95 ± 0.03 GeV obtained with the GKex(02S) model for 
the vector form factors from the fit to the subset of these data (160 datapoints). 
The data of Refs. [1,3,13-16,28] (39 grey datapoints) are rejected from the fit being 
either superseded or not satisfying our selection criterion x^/ndf < 4.5. 
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The resulting world average obtained are 

M^'^^[GKex(02S)] = 0.95 ± 0.03 GeV (xV^df = 0.92), (2) 

M^^'^[BBBA2006] = 0.96 ± 0.03 GeV (xV^df = 0.91). (3) 

The errors correspond to the usual one-standard-deviation errors (MINUIT 
default [42]) plus the systematic errors, added quadratically, which account 
for the uncertainties in the data on the vector form factors, nuclear effects 
(within the adopted model) and radiative corrections. The fit performed, 
for a comparison, with the naive dipole model for the vector form factors 
yields M^^^ = 0.93 ± 0.03 GeV with xV^df = 0.95. 

The obtained world average values 
of the axial mass are in strong con- 
tradiction with the recently published <«' 
result of the K2K Collaboration [44] : s 



M 



QES 



[K2K] = 1.20 ±0.12 GeV (4) 



This value has been determined for 

a water target through fitting the 

distributions of muon tracks 

reconstructed from neutrino-oxygen 

quasielastic interactions by using the 

combined K2K-I and K2K-IIa data 

from the Scintillating Fiber detector 

in the KEK accelerator to Kamioka 

muon neutrino beam.^ In Fig. |21 we 

cross section re- 
QES 




show the v^n 

calculated from the fit values of M 



K2K 2006 (recalculated from M^^ ) 
K2K / SKI default ( M*^= 1.10 GeV) 
K2K best-fit ( M?''^ = 1 .20 ± 0. 1 2 GeV) 
World average ( = 0.95 ± 0.03 GeV) 

'''''''' 3 ' ' 4 
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Fig. 2. Comparison between the QES 
^'/t^o — * cross sections (where no 
is a neutron bound in oxygen) evalu- 
ated with the best fit values JSJ and 10} , 
the current K2K/SKI default, and the 
K2K result reconstructed from the fit 
values of A/^^^ for five energy bins re- 
ported in Fig. 9 of Ref. [44]. 



obtained in Ref. [44] for the shape 
of the distribution for each recon- 
structed neutrino energy.^ The calcu- 
lation was performed with our default 
inputs that introduces an uncertainty 
of at most 2% which is added to the quoted error bars quadratically. Also 
shown are the cross sections evaluated by using the world average value © , 
the K2K best fit @, and the value of 1.1 GeV used as a default in the recent 
neutrino oscillation analyses of K2K [45,46] and Super-Kamiokande I [47]. 

^ Data from the continuation of the K2K-II period were not used in the analysis [44]. 

The best-fit values of M^^^ obtained from the K2K-I and K2K-IIa data subsets 

separately are, respectively, 1.12 ± 0.12 and 1.25 ± 0.18 GeV. 
^ The authors underline that the result for each energy should not be considered a 

measurement, but rather a consistency test. 



M'A printed on February 2, 2008 



5 



3. Axial mass from the data on single pion neutrinoproduction 

Figures I3I5I show a compilation of the data on single pion neutrinopro- 
duction cross sections from experiments at ANL [2,5], BNL [7], FNAL [8,9], 
CERN [18,19,21-26], and IHEP [31]. The nuclear targets are listed in the 
legends. All the data, as well as the theoretical curves, are classified through 
the panels corresponding to the experimental cut-offs in invariant hadronic 
mass W ranging from 1.4 to 2.55 GeV and including the measurements 
without cuts in W . 
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Fig. 3. Total tt and 7r+ production cross sections measured for different nuclear 
targets by the experiments ANL 1982 [5], BNL 1986 [7], GGM 1979 [21], BEBC 
1983 [23], BEBC 1990 [26], and IHEP SKAT 1989 [31]. The data are classified 
according to the cuts in W . The curves and bands correspond to the world average 
value of MY"^ = 1.12 ± 0.03 GeV obtained from the fit to a subset (196 points) of 
the full data presented in this and two next figures (see text for more details). 
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Fig. 4. Total tt , 7r°, and 7r+ production cross sections measured for different targets 
by tlie experiments ANL 1973 [2], ANL 1982 [5], BNL 1986 [7], FNAL 1978 [8], 
GGM 1978 [19], GGM 1979 [21], BEBC 1980 [22], BEBC 1986 [24], BEBC 1989 [25], 
BEBC 1990 [26], and IHEP SKAT 1989 [31]. See Fig. Hand text. 



M'A printed on February 2, 2008 



7 



• ANL 1982 (H,-DJ 



W<1.4GeV ; V^p ^[L^ nn" 

: * IHEP SKAT 1989 (CF,Br) 



W< 2.0 GeV 



V^n^H pn 

• ANL 1982 (H,-D,) 



W< 1.6 GeV 



V^^p^p. pT^ 

m GGM 1979 (CjHj-CF^Br) 
© BEBC 1983 (Dj) 
■ BEBC 1989 (Dj) 



W< 1.4 GeV 




V„n^\i pn" 



W< 2.0 GeV 



• ANL 1982 (H^-DJ, no cut 
BNL 1986 (Dj), no cut 

□ GGM 1978 (CjHj-CFjBr), W< 2.55 GeV 

* IHEP SKAT 1989 (CF^Br) 




V^p pn IV < 2.0 GeV 

m FNAL 1980 (H,), W< 1.9 GeV 

ffl GGM 1979 (CjHj-CF^Br), no cut 

A BEBC 1983 (Dj) 

A BEBC 1986 (Hj) 

* BEBC 1990 (Dj) 

* IHEP SKAT 1989 
(CF,Bi) 




10 



10 



10" 
(GeV) 



10 



10 

£„ (GeV) 



Fig. 5. Total tt" and tt production cross sections measured for different targets by 
the experiments ANL 1982 [5], BNL 1986 [7], FNAL 1980 [9], GGM 1978 [18], GGM 
1979 [21], BEBC 1983 [23], BEBC 1986 [24], BEBC 1989 [25], BEBC 1990 [26], 
and IHEP SKAT 1989 [31]. See Fig. Eland text. 



For the theoretical description of the single-pion neutrinoproduction 
through baryon resonances we apply an extended version of the Rein-Sehgal 
(RS) model [48]. Our extension [49] is based upon a covariant form of the 
charged leptonic current with definite lepton helicity and takes into ac- 
count the lepton mass. In the present calculations, we use the same set of 
18th nucleon resonances with central masses below 2 GeV and the same 
ansatz for the nonresonance background as in the original RS model. With 
that, all relevant parameters are updated according to the current data [41]. 
Significant factors (normalization coefficients etc.) estimated in Ref. [48] 
numerically are recalculated by using the new data and a more accurate 
integration algorithm. The relativistic quark model of Feynman, Kislinger, 
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and Ravndal [50] adopted in the RS approach unambiguously determines 
the structure of the transition amphtudes involved into the calculation and 
the only unknown structures are the vector and axial- vector transition form 
factors G^'^ (Q^). in the RS model, they are assumed to have the form 

with the "standard" value of the vector mass My = 0.84 GeV (that is the 
same as in the dipole parametrization of the elastic vector form factor). 
The integer n in the first ("ad hoc") factor of Eq. (jSJ is the number of 
oscillator quanta present in the final resonance. The axial mass Ma = 
M^^^ (fixed to be 0.95 GeV in the RS model) is the free parameter of 
our fit. In order to compensate for the difference between the experimental 
value of the nucleon axial-vector coupling qa and the SUq predicted value 
{gA{SU^) = —5/3), Rein and Sehgal introduced a renormalization factor 
Z = 0.75. For adjusting the renormalization to the current world averaged 
value gA = -1.2695 ± 0.0029 [41] we use Z = 0.762 and assume gy = 1. 

The nuclear effects for all nuclear targets different from hydrogen and 
deuterium are taken into account through the standard Pauli blocking fac- 
tor (see, e.g., Ref. [51] and references therein). The estimated relevant 
uncertainty is taken into account in the fit and in the error of its output. 

Almost all the data (196 points) shown in Figs. l3l5l DarticiDate in the fit. 
Several data subsets are excluded since they are superseded in the posterior 
reports of the same collaborations (e.g., the data from Refs. [2,23]), or are 
transformation of the others derived from the same experimental samples 
(e.g., the data of Refs. [21] with no cut on W). Note that all the data 
included into the fit satisfy the criterion x^/ndf < 4.5. The resulting world 
average obtained in the fit is 

M^^s = 1.12 ± 0.03 GeV (xV^df = 1.14). (6) 

As in the QES case, the error is the combination of the la deviation given 
by MINUIT and estimated systematic uncertainties. The obtained world 
average is in agreement with the recent analysis by Furuno et al. [52] of 
the BNL 7-foot bubble chamber deuterium data^ as well as with the value 

^ The analysis of Ref. [52] is based on the total event sample of 1.8 M pictures and 
holds two periods of runs in 1976-77 and 1979-80. The outputs of the analysis are 
MY'^ = 1.08 ±0.07 GeV (statistical error only) - from the fit of the distributions 
of pTT+ris events and — 1.15tQ;o6 GeV (both statistical and QES errors are 

included) - from the Itt and QES cross sections ratio. The best-fit value of 
obtained in the same analysis assuming the dipole model for the vector form factors 
(with the standard Mv) is 1.07±0.05 GeV that is well above our result (see Sect.|5J. 
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^RES = Qgv (the same as M^^^) adopted in the most recent K2K 
neutrino oscillation analysis [46] but considerably lower than the value 
= 1.2 GeV used for the atmospheric neutrino analysis of the Super- 
Kamiokande Collaboration [47]. 

Figure IHl shows a comparison between our calculations and the result 
of Ref. [52] for the ANL and BNL data on the ratios of the one-nucleon 
normalized Ivr and QES f^D2 cross sections (calculated and measured with 
no cut on W). Being transformations of the others, these data are not 
included into the fit. The narrow bands indicate the uncertainties in the 
values of the axial masses © and ©. The agreement is reasonably good. 
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Fig. 6. The ratios of Itt and QES cross sections evaluated in Ref. [52] from the data 
of ANL and BNL deuterium experiments. The curves and bands are calculated with 
the world average values 10) and ® for M^^^ and Mj^^^, respectively. 



4. Conclusions 

To summarise, we performed a statistical study of the QES and Ivr 
neutrinoproduction total cross section data in order to extract the best-fit 
values of the parameters M^^^ and . Our results given by Eqs. ©, 

(jSJ, and © are, of course, model dependent and can be recommended for 
use only within the same (or numerically equivalent) model assumptions as 
in the present analysis. We are planning to extend the analysis by employing 
more sophisticated treatments of the nuclear effects and including additional 
experimental information. 
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Appendix 

The most general formula for the QES i^N cross section is 

where s = {k + p)^, u = {k' — p)^, = —q^; k, k' = k — q, and p are the 
4-momenta of (anti)neutrino, final lepton, and initial nucleon, respectively; 
the coefficient functions A, B, and C are given by 

A=2[{x' + r2) (2.t' + _ Re (F^Fm) 

- {rRe [FX {Fy + Fm)] + [x + + x^) Re (F^Fp)} 
+ [{x' + (x' - 1 + - x^) - r^] [Fyp 

+ [{x' + (x' + 1 - - - r2] 
-[x' (x' + r2) + |FmP 

+ 4^2 (x' + x^) (x' + r^) |f|| 

± 4r (x + r^) [(x + 1 + x^) Re (F^Fa) + 2K^Re (F^Fp)] 
± 4rx2 [(x' + 1 + x^) Re (F|Fy) + x^Re (F|Fm)] 

- 4 (x' + r^) [(x' + x2) (x' + 1 + r^) + r^] |Fr|' 
+ 4x2 (^/^-L) (x' + x2) \Fs\\ 

B = ^ 4x'Re [Fl (Fy + Fm)] ± 2rx2 [|Fm|' + Re (F^Fa^ + 2FlFp) 

+ 4x2Re [F^ [Fa - 2 (x' + r^) Fp] - F^ {Fy - x' Fm) ] , 
C = |Fy p + |Fa|' + x' [Fm? T 4rRe {F^Fa) + 4 (x' + r^) |Ft|' , 

with the upper (lower) signs corresponding to neutrino (antincutrino) scat- 
tering. The six form factors Fj involved are functions of Q"^; 

,^ m ^ Mn - Mp ^ Mp + Mn 

^ 2{pq)' 4M2' 2M' 2M ' 2 ' 

and the remaining notation is standard. In the limit M„ = Mp, the general 

formula reduces to that of Ref. [40] and by putting F5 = Ft = it coincides 
with the result of Ref. [36] derived for the inverse (3 decay, taking account 
the proton-neutron mass difference.^ In this paper, we apply the Standard 
Model assumptions (T and C invariance + CVC), thus neglecting the scalar 
and tensor form factors Fs,t induced by the second-class currents, as well 
as the imaginary parts of the first-class form factors Fv,m,A,p- 



In fact, the latter efTect is insignificant for the present analysis, although it is included 
for completeness, together with the exact kinematics. 
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